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Experimental construction of a W-superposition state and its equivalence to the GHZ 

state under local filtration 
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We experimentally construct a novel three-qubit entangled W-superposition (WW) state on an 
NMR quantum information processor. We give a measurement-based filtration protocol for the 
invertible local operation (ILO) that converts the WW state to the GHZ state, using a register of 
three ancilla qubits. Further we implement an experimental protocol to reconstruct full information 
about the three-party WW state using only two-party reduced density matrices. An intriguing 
fact unearthed recently is that the WW state which is equivalent to the GHZ state under ILO, 
is in fact reconstructible from its two-party reduced density matrices, unlike the GHZ state. We 
hence demonstrate that although the WW state is interconvertible with the GHZ state, it stores 
entanglement very differently. 

PACS numbers: 03.67.Lx, 03.67.Bg, 03.67.Mn 


I. INTRODUCTION 

Explorations of multiqubit entanglement have un¬ 
earthed several families of states with curious quantum 
properties and there have been many attempts in recent 
years to characterize all the denizens of this quantum 
zoo [IH- The situation becomes complicated for sys¬ 
tems of more than two qubits and correspondingly the 
classification of their entanglement turns out to be more 
involved Bi- 

Pure entangled states of three qubits fall into two 
categories, namely the GHZ- or the W-class, under 
stochastic local operations and classical communication 
(SLOCC) [1, 0 with the maximally entangled GHZ and 
W states being given by: 

|GHZ) = ^(|000) + |111)) 

|W) = ^ (1001) + 1010)+ I100)) (1) 

The entanglement of the GHZ state is fragile under qubit 
loss, i.e when any one of the qubits is traced out^he other 
two qubits become completely disentangled [2, Isl . Hence 
if one of the parties decides not to cooperate, the entan¬ 
glement resources of the GHZ state cannot be used. In 
contradistinction to the GHZ state, the W-state residual 
bipartite entanglement is robust against qubit loss [2[. 

It has been shown by Linden et. ah, that almost ev¬ 
ery pure state of three qubits can be completely de¬ 
termined by its two-party reduced density matrices [9|. 
The two inequivalent entangled states, namely the W 
and GHZ states, have contrasting irreducibility features: 
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while GHZ states have irreducible correlations and can¬ 
not be determined from their two-party marginals 0 
ini, W-states are completely determined by their two- 
party marginals |I2l4l4j. Tripartite entan glem ent has 
been studied experimentally using optics llbl - H^ and 
NMR [Bil. 

Recently, the entanglement properties of a permuta¬ 
tion symmetric superposition of the W state and its ob¬ 
verse W = I/\/3(|100) + 1101) + 1110)) have been char¬ 
acterized [^, : 

|WW) = ^ (|W) + |w)) 

= 4 = (| 001 ) + | 010 ) + | 011 ) + | 100 ) + | 101 ) + | 110 )) 
v6 

( 2 ) 

While this state (referred to henceforth as the WW state) 
belongs to the GHZ entanglement class, its correlation in¬ 
formation (in contrast to the GHZ state) is uniquely con¬ 
tained in its two-party reduced states. The argument for 
reconstructing the three-qubit WW state from its two- 
party reduced states runs along similar lines to the orig¬ 
inal argument of Linden et. al. [9|. If we assume another 
state to have the same two-party reduced density matri¬ 
ces as the WW state, this constraint can be used to prove 
that the new state is no different from the original WW 
state m, [26|. 

In this work we focus on the WW state. We provide an 
explicit measurement-based filtration scheme to filter out 
the I GHZ) state from the WW state. Further, we exper¬ 
imentally construct and tomograph the WW state on an 
NMR quantum information processor of three coupled 
qubits. We experimentally demonstrate that the infor¬ 
mation about tripartite correlations present in this state 
can indeed be completely captured by its two-party re¬ 
duced density matrices. We reconstruct the experimental 
density matrices using complete state tomography and 
compare them with the theoretically expected states and 
also compute state fidelities. The GHZ class of states are 
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an important computational resource [l| and it has been 
shown that states that are SLOCC equivalent to these 
can be used for the same kind of quantum information 
processing tasks [2|. Therefore, it is expected that the 
WW state will also prove useful for quantum computa¬ 
tion. Furthermore, the quantification of the tripartite 
correlation information present in this state is easier as 
compared to the GHZ state, as entanglement measure¬ 
ment requires only two-qubit detectors. 

The paper is organized as follows: Section HIl describes 
how we obtain the GHZ state from the WW state by lo¬ 
cal filtration based on projective measurements using a 
register of three ancilla qubits. Section [1111 describes the 
experimental creation of the WW superposition state on 
a three-qubit NMR quantum information processor. Sec¬ 
tion IIII Al contains the details of the molecule used, the 
NMR pulse sequence for WW state construction and the 
results of state tomography. The information content of 
the WW as captured from its two-party marginals is de¬ 
scribed in Section illlBl We conclude in Section ilVl with 
some remarks about GHZ and WW types of three-qubit 
entanglement and the relationship between entanglement 
class and how information about entanglement is stored 
in a quantum state. 


II. FILTRATION PROTOCOL TO SHOW 
SLOCC EQUIVALENCE OF WW AND GHZ 


Measurement-based local filters have been used for en¬ 
tanglement manipulation in the context of violation of 
Bell inequalities as well as for the detection of bound en¬ 
tangled states No local operations can convert 

a state from the GHZ class to the W class. However, 
surprisingly, it has been shown that the WW is in the 
GHZ class, deriving from the fact that it is related to 
the GHZ state via the SLOGG class of operations given 
by @,111: 


\GHZ) = A (8) A (8) A IVFVF) (3) 

with 


A 




(4) 


being an ILO, where uj = denotes the cube root of 
unity. We have used ‘=’ instead of an equality sign in 
Eqn. ([3]) because A is a non-unitary operator that does 
not preserve the norm and the two sides in Eqn. ([3]) do 
not have the same norm. 

We now proceed to reinterpret A as an action on an 
ensemble of identically prepared states WW and imple¬ 
ment the operation described in Eqn. ([3]). In this pro¬ 
cess, we will have to discard some copies and the new 
ensemble that we construct with each member in the 
filtered GHZ state will have fewer copies as compared 
to the original ensemble of WW states. These aspects 
will be brought out more clearly when we describe the 
measurement-based filtration protocol to realize the ILO. 


Since A acts on each of the qubits locally, we first want 
to realize the operation A on a single qubit. The non- 
unitary operator A has a singular valued decomposition 


A = UDV (5) 


where the unitary operators U and V are given by 






(6) 


and the non-unitary diagonal operator D is given by 


-Ci) 

The operators U and V are unitary and can be imple¬ 
mented via a local Hamiltonian evolution. Therefore, 
we now turn to the implementation of D on a one-qubit 
state. 

Erom the two columns of the operator D we define two 
vectors 


K)=(J) and M = ^(3") (8) 


These vectors are orthogonal to each other but are not 
normalized. We now extend the Hilbert space of the 
system by adding an ancilla qubit. We extend the vectors 
ui and U 2 to the composite Hilbert space formed by the 
ancilla and the system to obtain two four-dimensional 
vectors 


iei) = 


and IC 2 ) = 


^/3 


0 

33 

0 

Vys - v^j 


(9) 


The vectors |^i) and 1 ^ 2 ) are not only mutually orthogo¬ 
nal but also normalized. 

Using these orthonormal vectors |^i) and 1 ^ 2 ), we con¬ 
struct orthogonal projectors Pi and P 2 


i"i = iy)(yi 

P 2 = IC 2 ) (61 


/I 0 0 0\ 
0000 ] 
0 0 0 0 
Vo 0 0 0/ 


1 

7 ! 


/o 0 0 0 \ 

0 1 0 yvs -1 

0 0 _ 0 0 

Vo yys-10 vs-i / 


( 10 ) 


We define the projection operator P = Pi -f P 2 - The 
effect of the projector P on the composite system of the 
single qubit and a one-qubit ancilla turns out to be 


P = 


' D 

A ' 

^ A 

D' 


4x4 


(11) 
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FIG. 1. Schematic diagram of the filtration scheme to im¬ 
plement the non-unitary ILO transformation that converts a 
WW state to a GHZ state. 


where D is the diagonal part of the singular value decom¬ 
position of the operator A given in Eqn. m, the comple¬ 
mentary matrix D' = I — D and the matrix A can be 
obtained readily from Eqn. ([TO]) . 

If we prepare the ancilla in a state |0) (0| with the 
system being in an arbitrary state p, the action of P 
on the composite system is given by 


P(|0) (0|(8)p)P 


DpD 

Dp A ' 

ApD 

<1 

<1 


( 12 ) 


If we measure the projector P on the composite system 
(system and ancilla), whenever the measurement gives 
a positive answer, the state after measurement is given 
by the right hand side of Eqn. m- We retain only 
these cases and discard the state whenever the outcome 
of the measurement is negative. Eurther, on the final 
state given in Eqn. m, we measure the projector |0) (0| 
on the ancilla alone. As before, if the outcome is positive 
we retain the state, and if the outcome is negative we 
discard the state. In case the outcome is positive, the 
resultant state is |0) (0| (8) DpD and upon discarding the 
ancilla we get the state of the system to be DpD. This 
completes the application of the non-unitary invertible 
operator D on p. Sandwiching this operation between the 
unitary transformations U and V as given in Eqn. 
we achieve the application of the ILO operator A on p. 

The scheme is easily extendable to 2 O 2 (8) 2 systems, 
where we locally implement A on each of the three qubits. 


We imagine that the tripartite system is divided between 
Alice, Bob and Charlie and each of them can perform 
local operations at their location. We begin with the 
state |WW) for the three qubits, attach a one-qubit an¬ 
cilla to each qubit, and measure the local projector P for 
each qubit. If the outcome of these measurements (that 
amount to a measurement of P (8) P G) P) is positive we 
retain the state, otherwise we discard the state. Then on 
each ancilla, we measure the projector |0) (0| and retain 
the cases when all the outcomes are positive. Upon dis¬ 
carding the ancillas, the resultant state is the application 
of D on each qubit. When we sandwich this process be¬ 
tween the unitaries U and V on each qubit, we get the 
final state as |GHZ). This process of measurement-based 
filtration is schematically explained in Eig. [H To decide 
when to discard and when to retain the outcome, we re¬ 
quire classical communication between Alice, Bob and 
Charlie. Since we discard the output state in a number 
of cases, the size of the ensemble obtained in the end is 
smaller than the original ensemble. 


III. NMR IMPLEMENTATION 



=9497.0 Hz 
1^2=0.0 Hz 
iy3=-13897.0Hz 


Ji2=69.8Hz 

J23=-129.0Hz 

Ji3=47.5Hz 


(b) 


Spin 1 Spin 2 Spin 3 

111 ) | 10 ) | 01 ) | 00 ) | 10 ) | 00 ) 111 ) | 01 ) | 10 )| 00 ) | 11 )| 01 ) 



(jJf (in ppm) 


FIG. 2. (a) Molecular structure and NMR parameters 

(chemical shifts and J-coupling in Hz) and ^^F NMR spectrum 
of trifluoroiodoethylene. The three fluorine spins correspond 
to the three-qubit system, (b) The ID ^®F NMR thermal 
equilibrium spectrum obtained after a ^ readout pulse. The 
NMR transitions of each qubit are labeled by the correspond¬ 
ing logical states of the other two qubits. 


To prepare the WW state on a three-qubit NMR quan¬ 
tum information processor, we employ the three fluorine 
(spin-1/2) qubits of trifluoroiodoethylene. The molecular 
structure and NMR parameters of this three-qubit sys¬ 
tem are adequate for the kind of manipulations involved 
in quantum state preparation and are given in Fig. [2fa). 
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Average fluorine longitudinal Ti relaxation times of 5.0 s 
and T 2 relaxation times of 1.0 s were experimentally de¬ 
termined. The equilibrium fluorine NMR spectrum ob¬ 
tained after a ^ readout pulse is shown in Fig. EJb). 

The system was first initialized into the |000) pseu- 
dopure state using the standard spatial averaging tech¬ 
nique [ 3 ^. The experimental density matrices were 
tomographed by standard state tomography proce¬ 
dures [3ll-[^ . The three-qubit experimental density ma¬ 
trix was tomographed using a set of eleven detection op¬ 
erators defined by {III, IIX, IXI, XII, IIY, lYI, YII, YYI, 
IXX, XXX, YYY}, and the two, two-qubit reduced den¬ 
sity matrices were determined using a set of four detec¬ 
tion operators defined by {III, IXI, lYI, XXI} and {III, 
IIX, IIY, IXX} respectively, with I denoting the identity 
(or no-operation) operator and X(Y) denoting a spin- 
selective ^ pulse of X(Y) phase on a specified qubit. The 
fidelity of the reconstructed state was computed using 
the Uhlmann-Jozsa fidelity measure 00: 


F = 


(^'sj \/PtheoryPexpt\/Ptheory^ ^ 


where ptheory and Pexpt denote the theoretical and exper¬ 
imental density matrices respectively. 


A. WW construction scheme 


The circuit to construct a WW state consists of sev¬ 
eral single-qubit and two-qubit gates. A single-qubit gate 
Fi[a]y acting on the ith qubit, achieves a rotation by the 
angle a around the y axis with a corresponding unitary 
matrix given by: 


Ui[(y\y — 


cos ^ — sin ^ 

• cy. cy 

sm ^ cos ^ 


(14) 


A two-qubit controlled-rotation gate CRij[0]y, imple¬ 
ments the single-qubit rotation Uj[(l)]y on the target qubit 
j about the y axis, if the control qubit i is in the state 
11). The CNOTij gate implements a controlled-NOT op¬ 
eration with the ith qubit as control and the jth qubit 
as target. 

The sequence of gates to construct a WW state, start¬ 


ing from the initial pseudopure state |000) is given as: 


| 000 ) 





1 




i (\/3|000) - |100)) 




CR 12 

[2 cos 

1 


i (\/3|000) - 


0 ^ 

1 

1 







i (\/3|000) 

-73 (!“»> + 

■ |110) + 

1010))) 



CNOT 13 




i (y3|000) 

-^(1101) + 

■ 1111 : 

) + 

1010))) 



CNOT 23 




i (V3|000) 

-^(|101) + |110) + 

1011))) 

1 

Ui[f 


3[f]. 

1 



^ (| 001 ) + 1010 ) + 1011 ) + | 100 ) + | 101 ) + | 110 )) 

The quantum circuit to construct the WW state on a 
three-qubit system is given in Fig.[3}a). The NMR pulse 
sequence to create the WW state, starting from the pseu¬ 
dopure state 1000) is given in Fig. [3}b). All the pulses 
are shaped pulses, labeled by the corresponding axes of 
rotation and the flip angles; Tij denotes an evolution 
period under the Jij coupling. Refocusing (tt) pulses 
are applied in the middle of the evolution periods to 
compensate for chemical shift evolution and pairs of tt 
pulses are introduced at 1/4 and 3/4 of the evolution 
periods to eliminate undesired J-evolutions. After the 
evolution interval r 23 and the [^]y on the third qubit 
(corresponding to a CNOT 23 gate), the state obtained is 
^|000) — + |110) + ^|0II)). There is an un¬ 

desirable extra relative phase of h’ that has accumulated 
between two of the basis vectors. This undesired extra 
phase factor is compensated for during the evolution in¬ 
terval 3 ri 2 . The implementation of the last module (si¬ 
multaneous [^]y pulses on all the three qubits) results in 
the desired WW state with no extra relative phase. All 
the selective pulses are 265 /is Gauss’^ shaped pulses 
and the non-selective excitation pulse is a frequency- 
modulated 400 /is ^^Gauss^^ shaped pulse. The NMR 
spectrum of the WW state obtained by a sequence of se¬ 
lective rotations on the initial pseudopure state is shown 
in Fig. m Each spin multiplet has two resonance peaks 
(as compared to four resonance peaks for the thermal 
equilibrium state). The expected NMR spectral pattern 
of an ideal WW state should contain resonance peaks 
of equal magnitude and phase, and deviations from ideal 
spectral peak intensities and phases in the experimentally 
obtained spectrum, can be attributed to imperfections in 
the rf pulse calibrations and to relaxation during the se¬ 
lective pulse durations. 
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(a) 



FIG. 3. (a) Quantum circuit showing sequence of gates re¬ 
quired to construct the WW state, starting from the pseu- 
dopure state |000). The gate operations are described in 
the main text and all the rotations take place about the 
y-axis. (b) NMR pulse sequence to create a WW state. 
All the pulses are low-power selective pulses represented by 
shaped blocks. Filled black shapes are tt refocusing pulses, 
unfilled shapes correspond to pulses of ^ flip angle and the 
gray shaded shapes are labeled with their specific flip an¬ 
gles and phases. The axes of rotation are specihed at the 
top of each pulse. Vertical dotted red lines show the corre¬ 
spondence between the quantum circuit and the experimen¬ 
tal pulse sequence. All pulses are of phase x unless other¬ 
wise labeled. The values of the rf pulse flip angles used are 
a — ^ — 2 cos “^(^),7 = ^ and rij represents an evo¬ 

lution under the Jij coupling. The last 3 ti 2 period is used 
to compensate the extra phase acquired (as described in the 
text). 


The tomograph of the experimentally constructed WW 
state is shown in Fig. O The experimentally tomo- 
graphed state was compared with the theoretically ex¬ 
pected state and the density matrices match well, within 
experimental error, with a computed state fidelity of 0.94 
(the fidelity was computed from Eqn. [T3|) . 


B. Reconstruction of WW from two-party reduced 
density matrices 

A protocol was developed [T2| to validate the surpris¬ 
ing aspect of multi-party correlations asserted by Linden 
et. al. [9|,[36|, that the information about three-party cor¬ 
relations of almost all pure three-qubit states (except for 


Spin 1 Spin 2 Spin 3 

111 ) | 10 ) | 01 ) | 00 ) | 10 ) | 00 ) 111 ) | 01 ) | 10 )| 00 ) | 11 )| 01 ) 



FIG. 4. The ID ^®F NMR spectrum corresponding to the 
creation of the WW state. The NMR transitions of each qubit 
are labeled by the corresponding logical states of the other two 
qubits. 


Re Im 





FIG. 5. The real (Re) and imaginary (Im) parts of the (a) 
theoretically expected and (b) experimentally density matri¬ 
ces for the WW state reconstructed using full state tomog¬ 
raphy. The rows and columns of the bar graphs depict the 
computational basis of the three qubits in binary order from 
|000) to I 111). The experimentally tomographed state has a 
fidelity of 0.94. 


GHZ-type states) are already contained in their corre¬ 
sponding two-party reduced states. We delineate below 
the argument for how a general three-qubit pure state 
pABC can be completely determined by using any of the 
equivalent sets {pab,Pac), {pab,Pbc), or {pac,Pbc) 
of reduced two-party states. The reduced single-qubit 
reduced state pA and the two-qubit reduced state pBC 
share the same set of eigen values, and can hence be 
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written as 


Pa = Y1pa \i) {i\ 

i 

PBC = Y.PA\i-,BC) {v,BC\ 


(15) 


where {|i)} are the eigenvectors of pA with eigenvalues 

eigenvectors of pBC with 
eigenvalues {p^}. Furthermore, the three-qubit pure 
states that are compatible with pA and pBC are given 
by: 


\'tpA 


BC] <a 


) = E 


|i) (8) \i; BC) (16) 


Similarly, the three-qubit pure states that are compatible 
with pc and pab are given by 

\'>PABC-n) = Y.^^''"\[^c\k\AB)®\k) (17) 

k 

where {\k)} are the eigenvectors of pc with eigenval¬ 
ues {pj} and {\k]AB)} are the corresponding eigen¬ 
vectors of paB‘ Since the pure state is com¬ 

patible with both Pab and pBCi we can now consis¬ 
tently find the values of ai and y/c while ensuring that 
|'0ABC;<a) = \i^ABC]l)- 

We used the set of two, two-party reduced states 
{pABi pBc)i to reconstruct the full three-qubit WW 
state. The reconstructed density matrix for the WW 
state, using two sets of the corresponding two-qubit re¬ 
duced density matrices (pab^ Pbc) is given in Fig. [6l 
The two-party reduced states were able to reconstruct 
the three-party WW state with a fidelity of 0.92, which 
matches well with the full reconstruction of the entire 
three-qubit state given in Fig. [5{b). 


IV. CONCLUSIONS 

We described a measurement-based filtration scheme 
to demonstrate the ILO equivalence of the WW state 
with the GHZ state. We experimentally implemented 
an NMR-based scheme to construct a WW state. We 
were able to show that the three-qubit density operator 
PABC obtained by full state tomography matches well 
with the same three-qubit state reconstructed using a 
set of two-party reduced density operators {pab^ Pbc)- 
Thus, although the WW state belongs to the same en¬ 
tanglement class as the GHZ state, the two states store 
information about multi-party correlations in completely 
different ways. We thus experimentally demonstrated an 
interesting feature of multi-qubit entanglement namely, 
that two different entangled states belonging to the same 
SLOGG class can yet have their correlations exhibiting 
contrasting irreducible properties. 

Since distinguishing entangled states is still a hard 
task, our work can be used as a benchmark to further 


(a) Pab 

Re Im 


(b) PBC 

Re Im 


® 


Ot xy 


(c) 


PABC 


Re 


Im 



FIG. 6. The real (Re) and imaginary (Im) parts of the to- 
mographed density matrix for the state WW state, (a) The 
two-qubit reduced density matrix pab • (b) The two-qubit re¬ 
duced density matrix pbc- ( c ) The entire three-qubit density 
matrix pabc, reconstructed from the corresponding two-qubit 
reduced density matrices. The rows and columns in the bar 
graphs encode the computational basis of the qubits, from 
100) to 111) for two qubits and from |000) to |111) for three 
qubits. The fidelity between the three-qubit state (pABc) 
reconstructed from the two-qubit density matrices and the 
three-qubit state obtained by complete three-qubit state to¬ 
mography is found to be 0.92. 


classify how different entangled states store information 
about their correlations. Our work also has important 
implications for comparing the utility of different kinds 
of entangled states to perform the same computational 
task. We were unable to find a suitable molecular ar¬ 
chitecture to experimentally implement the ILO, since 
this requires each of the three qubits to be coupled to a 
separate one-qubit ancilla. However, it is a worthwhile 
exercise to look for an experimental implementation of 
the filtering protocol to perform the ILO. A further issue 
with such an implementation is the involvement of pro¬ 
jective measurements, which are not straightforward to 
achieve using NMR. 
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